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THE  CANTED  SPECTRAL  TRANSFORM 
AND  ITS  PROPERTIES 


INTRODUCTION 


Spectral  transforms  play  a  major  role  in  modern  signal  processing  and  analysis.  This  role  will 
undoubtedly  be  expanded  in  the  future  as  a  consequence  of  advances  in  microelectronics  and  computer 
technologies  The  most  common  spectral  transform  used  in  practice  is  the  discrete  sectionalized 
Founer  transform  (SFT)  which  evolved  with  the  advent  of  the  fast  Fourier  transform  (FFT)  algorithm 
of  Cooley  and  Tukey  [1,2].  The  FFT  algorithm  opened  the  door  to  high-speed  coherence  estimation  in 
ow-frequency  (acoustic)  applications  through  its  use  as  an  SFT  [3-6],  In  application,  the  SFT  accumu- 
ates  a  signal  energy  along  narrow  spectrally  invariant  channels  (or  frequency  bins)  over  the  temporal 
hmits  of  the  transform  integration  interval.  This  restriction  on  signal-spectral  dynamics  is  imposed  by 
the  Fourier  kernel,  exp  {-i2nft}.  One  may  therefore  surmise  that  this  restriction  is  rather  arbitrary 
and  that  one  may  accumulate  signal  energy  along  any  dynamic  spectral  path  by  choosing  a  suitable  ker¬ 
nel  function  to  use  in  the  transform.  Suppose,  for  example,  that  one  chooses  a  kernel  exp  {-/¥(*)} 
where 


ITT  2^  Vn*n  +  l  “  27T  f  7)0  t  +  T)\  t2  +  .  .  .  4*  tN  1  . 

*=o  1  J 


(la) 


The  transform  of  a  signal  using  W(t)  as  a  kernel  will  accumulate  energy  over  dynamic  spectral  paths 

/<>>-■£- V  (t)  -  X>+1)„./".  (lb) 

12— 0 


Here  'F(r)  is  the  time  derivative  of  the  generalized  phase  function  V(t).  It  should  be  evident  there¬ 
fore,  that  for  a  narrowband  signal  s(()  -  A  sin  0  it)  whose  instantaneous  frequency  satisfies 


1  ’  /  x  N~l  (0) 

9(t)  =  T  J- - — 

llT  n- 0  »! 


(2a) 


over  0  <  t  <  T,  the  magnitude  square  of  its  spectral  transform  reduces  to 


1_ 

T 


=  \  So  s2  {t)dt . 


(2b) 


for 


rin  =/(n)  (0)/  («  +  !)! 


(2c) 


With  foreknowledge  of  the  spectral  dynamics  of  a  given  signal,  the  kernel  parameters  Vn  of  the  general¬ 
ized  transform  may  be  chosen  to  provide  an  effective  matched  filter  for  the  signal.  Without  this  fore- 
nowledge,  the  kernel  parameters  may  be  varied  to  search  for  the  set  that  maximizes  the  'transform 
output;  thus  providing  an  estimate  of  the  detected  signal  spectral  dynamics. 
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Although  sound  in  principle,  the  generalized  spectral  transform  becomes  increasingly  computa¬ 
tionally  intensive  as  the  order  of  the  phase  polynomial  becomes  large.  However,  extending  the  phase 
kernel  from  linear  to  quadratic  is  a  reasonable  step  to  further  the  cause  of  the  generalized  concept  in 
spectral  transforms. 

The  utility  of  the  quadratic  phase  kernel  has  recently  received  some  degree  of  recognition.  Wol- 
cin  [7]  has  demonstrated  that  when  the  spectral  trajectory  of  a  signal  can  be  adequately  .modeled  by  a 
continuous,  piecewise  linear  function  of  time,  the  maximum  a  posteriori  estimate  of  the  narrowband 
signal  (not  surprisingly)  results  in  exponential  quadratic  phase  functions.  George  Rogers,  at  TRW, 
Inc.,  subsequently  determined  that  the  use  of  the  quadratic  phase  kernel  in  place  of  the  conventional 
linear  Fourier  kernel  could  significantly  improve  the  detection  of  particular  classes  of  spectral  dynamic 
signals  propagating  in  the  ocean  medium  [8].  The  objective  of  this  report  is  to  study  the  properties  of 
the  spectral  transform  that  employs  the  quadratic  phase  kernel,  in  a  preliminary  effort  to  determine  its 
utility  in  practical  applications.  This  form  of  the  spectral  transform  will  be  called  the  canted  spectral 
transform  (CST);  in  that  it  accumulates  signal  energy  along  spectral  paths  that  may  be  canted  (or 
linearly  sloping)  with  respect  to  time.  The  discrete  form  of  the  transform  using  the  combined  Fourier- 
Fresnel  kernel  is  also  referred  to  (by  Rogers)  as  the  slide-fast  Fourier  transform  or  S-FFT. 

FORMULATION  OF  THE  CST 


The  discrete  form  of  the  CST  to  be  studied  in  this  report  is* 


F{k%v) 


\_ 

W 


N- 1 


—  i2n 


L  Wn\+n  e 


N 


2  N 


(3) 


where 

Sj  are  the  signal  samples  over  the  interval  J  <  ~F  TV  —  1 , 
k  =  0,1, 2, 3, . . .  ,N—  1  is  the  frequency  index, 

z^  =  0,±l,  ±2,  ±3,  . . .  ,±z/max  is  the  cant  (or  frequency-slide)  index, 
n  is  the  time  index  (or  sample  number)  relative  to  «0, 
n0  is  the  initial  time  index, 

wn  =  wN_n  is  the  window  function  over  0  <  n  <  N-\, 


N- 1 

W  =  £  w„  is  the  windowed  signal-averaging  factor. 

n= 0 

The  signal  s(t )  is  assumed  to  be  uniformly  sampled  over  time  increments  A/.  The  transform- 
integration  time  T  is  therefore  N A  t.  The  bin  width  of  the  spectral  channels  is  1/2T,  centered  at 
/  =k/T.  The  frequency  slide,  f  —  A//  T  =  A/c/  T2,  of  the  spectral  channels  is  v/  T 2.  The  cant  v  is  the 
integer  number  A  A;  of  frequency  bins  shifted  over  the  transform-integration  time  T.  The  factor  W  is 
equal  to  N  for  the  conventional  rectangular  window  and  equals  N/2  for  the  Hanning  window.  The 
spectral  channel  width  1/  T  is  generally  chosen  to  encompass  the  short-term  bandwidth  of  the  nar¬ 
rowband  signal  components  of  interest.  And  the  range  of  the  cant  variable  is  chosen  to  encompass  the 
anticipated  frequency  shift  of  the  signal  components  over  the  transform-integration  time.  As  a  conse¬ 
quence,  the  class  of  narrowband  signals  most  suitable  for  analysis  by  the  CST  is  that  whose  member’s 
instantaneous  frequency  deviates  from  linear  over  time  T  no  greater  than  approximately  1/  T. 


*Although  the  discrete  CST  would  be  complex  and  not  include  the  absolute  value  signs,  it  is  convenient  to  consider  only  the 
magnitude  of  the  transform  in  this  early  study  of  its  general  properties. 
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THE  CHARACTERISTIC  FUNCTION  Z(x,y) 

Since  the  CST  is  structured  to  accumulate  signal  energy  along  canted  or  sloping  channels  in  the 
frequency-time  plane,  let  s(t )  be  a  signal  whose  instantaneous  frequency  varies  linearly  with  time. 
That  is,  let 


=  e 


27rk0-jv+'trv0 


N 


+  Or 


(4) 


where 

*o  =  foT 

fo  is  the  instantaneous  frequency  of  s(t)  at  t0  =  noAt 
vo  =  foT2 

/o  is  the  time  derivative  of  the  instantaneous  frequency  at  t$ 


and 


0q  is  a  constant  phase  term. 


In  this  event,  from  Eqs.  (3)  and  (4),  the  CST  becomes 


F (kyp)  = 


w 


N- 1  -iir 

YW»e 

rt-0 


Letting  x  =  k—  kQ  and  y  =  v  —  defines  the  characteristic  function  of  the  CST  as 

)2] 

n  m 

N—\ 

Y,  wn  e 


1 


ZU,J,)  =  1 V 


2^+- 


N 


n-  0 


(5) 


(6) 


The  CST  characteristic  function  Z(x,y)  is  a  quasi-continuous  function  of  the  variables  x  and  y  (since 
k0  and  p0  need  not  be  integers),  and  it  serves  to  demonstrate  some  of  the  properties  of  the  CST. 

Symmetry  Relations 

Replacing  x  by  -x  and  y  by  -y  in  Eq.  (6)  gives  the  symmetrical  relationship 

Z(-x,-y)  =  Z(x,y) .  (7a) 

The  above  symmetry  implies  that  paired  values  k,  k'  and  v,  v'  (corresponding  to  pairs  x,-x  and  y,-y 
respectively)  must  satisfy  the  relations 

k  =  (k  +  k')/2  =  k0,  (7b) 


and 


v  =  (v  +  v')/2  =  vo  •  (7c) 

Consequently,  k0  and  v0  must  be  integer  multiples  of  1/2  since  k,  k\  v ,  and  v' are  integers.  The 
geometric  interpretation  is  that  k  and  k '  are  spectral  bins  symmetrically  located  about  k0  with  cants  v 
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and  v '  such  that  the  resulting  spectral  trajectories  intersect  the  spectral  trajectory  of  the  signal  k0  vn  at  a 
common  point. 


Note  from  Eq.  (5)  that  the  CST  is  dependent  on  only  the  differences  k  -  k0  and  v  - 1>0.  Thus  one 
can  determine  the  canted  transform  of  any  linear  frequency-slide  signal  (including  the  cw  signal)  by 
computing  the  conventional  Fourier  transform  of  a  canted  signal  whose  parameters  are  k0  and  v0-v. 
This  is  apparent  from  the  form  of  the  CST  in  Eq.  (3),  which  can  be  interpreted  as  the  conventional 
Fourier  transform  of  a  signal  modified  by  the  phase  kernel,  exp  {-iirvn2/ N2}. 


Letting  n  go  to  N  n  in  Eq.  (6),  and  using  the  identity  =  w„,  the  CST  characteristic  func¬ 
tion  becomes 


Z(x,y ) 


W 


N- 1  hr 


£ 

n— 0 


2 

nx+y)jf-y 

n 

N 

H>0 


]  __  ein(2x+y) 


(8a) 


For  typically  large  values  of  N ,  the  second  term  is  insignificant  in  comparison  to  the  summed  term  and 
can  be  ignored.  Therefore, 

Z(—x—y,y)  =  Z(x,y)  -  Z(x+y-y)  =  Z(-x,-y)  .  (8b) 

Thus,  the  two  values  k  and  k '  =  2k0  +  v0—  v  —  k  (corresponding  to  x  and  —x~y)  will  give  the  same 
value  for  the  characteristic  function.  In  geometric  terms,  given  the  values  k0,  v0,  and  v,  any  two 
values  of  k  whose  arithmetic  mean  is 

k  =  (k  +  k')/ 2  =  k0+  (vq—  v)/2  (8c) 


will  give  the  same  value  for  the  characteristic  function.  Writing  the  above  equation  as 

k  +  v/2  =  ko  + vq/2  (8d) 


reveals  that  the  spectral  trajectory  originating  at  k  with  cant  v  intersects  the  trajectory  originating  at  Ar0 
with  cant  v0  at  its  midpoint.  Consequently,  the  trajectories  originating  at  k  and  k '  will  intersect  the  tra¬ 
jectory  of  the  relevant  signal  s  (r)  at  points  symmetrically  distributed  about  its  midpoint.  Since  k  and  v 
are  integers,  it  is  necessary  that  2k0  +  v0  be  an  integer. 

Translation  of  the  Temporal  Index 


Translating  the  temporal  index  of  the  CST  characteristic  function  by  a  gives  the  interesting  result 

that 


Z'(x,y;a)  =  jp 


W 


N- 1  -hr 

£  e 

n=Q 


N~  1  -I* 

£  e 

n=Q 


2x 


+  y 


N 


(9a) 


=  Z ( x-ay,y ) . 


Thus,  translating  the  time  index  n  by  a N  is  equivalent  to  translating  the  x  variable  by  ay.  In  particu¬ 
lar, 
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Z'(0,y-x/y)  =  Z(x,y) 


_1_ 

W 


T 

T 

_n_  2£_ 

2 

L  wn  e 

N  y 

n= 0 

(9b) 


Thus,  one  can  compute  the  characteristic  function  using  only  the  quadratic  phase  term  by  appropriate 
translation  of  the  time  index. 


A  study  of  Eq.  (9b)  reveals  that  the  range  of  phase  excursion  over  n  is  minimized  when 
x  =  —y/2.  The  characteristic  function  for  minimum  phase  excursion  may  be  written  as 

>2 


Z(-y/2,y)  =  ± 


2_ 

W 


N/2 

2  j  WN/2+n  e 
n-0 


-  [  1  +  w0e~ll,y/* j 


2 

N/2  -iny 

n 

N 

X  WN/2  +  n  e 

n-0 

(9c) 


Topology  of  Z(x,y) 


To  map  the  topology  of  the  CST  characteristic  function  over  the  x,  y  plane,  the  functional  form 
given  in  Eq.  (9b)  may  be  used  to  give 


where 


Z(x,y)  —  |  ZR  (x,j/)  +  Zf(x,y)  j  , 


1  N—  l 

ZR(x,y)  =  — }  Y,  wn  cos 


n= 0 


rry 


n  x 
N  y 


c 1 

2 

J0  >P(f)  COS 

Try 

y 

d(9 


(10a) 


(10b) 


and 


zi(x,y)  -  -jp  £  w„  sin 


n= 0 


Try 


n  x 
N  y 


f  w(f)  sin 

iry 

e+- 

2' 

J  u 

y  : 

dZ. 


(10c) 


Although  the  above  functions  are  undefined  for  y— 0,  the  singularity  is  removable,  resulting  in  the 
form  of  the  conventional  Fourier  transform  of  the  window  function. 

Rectangular  Window 

For  a  rectangular  window  function,  a  closed  form  of  the  solution  is  obtained  involving  the  Fresnel 
sine  and  cosine  integrals  S  ( •  •  • )  and  C  ( •  •  • )  as  follows: 


$2 


(11a) 
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Z,{x,y)  = 


S  (p2  +  p ,)-$  (ffi) 

Pi 


where 


P\  =  ± 


Ul  and  p2  =  V2T7T, 


(lib) 


(11c) 


and  where  the  sign  of  /Sj  is  positive  when  x  and  y  have  the  same  sign  and  negative  when  their  signs 
differ.  This  result  for  a  linear  slide  frequency  was  demonstrated  in  Ref.  6  and  has  also  been  derived 
and  demonstrated  by  Rogers  [8]  in  his  unpublished  work  on  the  slide-FFT.  For  x=0  or  x=-y,  the 
characteristic  function  reduces  to 


Z(0,y)  =  Z(-y,y)  = 


v  s2(fi2) + c2(/j2) 
Wi 


(12a) 


and  for  x=—y/2. 


Z(-y/2,y) 


Vs2(/32/2)  +  C2(j32/2) 

Pi/2 


(12b) 


Figure  1  shows  graphical  plots  of  the  CST  characteristic  function  for  the  rectangular  window.  To 
achieve  the  realism  of  the  CST,  the  topology  is  mapped  over  integer  increments  of  x  and  y  about  initial 
offsets  of  either  0  or  1/2.  As  a  visual  aid,  the  discrete  values  along  the  x  axis  are  connected  by  straight 
lines.  The  offsets  of  1/2  in  either  or  both  the  x  and  y  axis  are  intended  to  depict  the  effect  of  nonin¬ 
teger  values  for  the  signal  parameters  k0  and  p0  on  the  resulting  topology.  In  the  upper  left-hand 
diagram,  the  parameters  k0  and  p0  are  both  integers.  Consequently,  a  peak  value  of  1  is  achieved  for 
x=y= 0.  In  the  upper  right-hand  diagram,  the  parameter  k0  is  depicted  to  be  an  odd  multiple  of  1/2. 
That  is,  the  initial  frequency  is  depicted  to  start  at  the  edge  of  a  CST  bin,  while  the  cant  p0  remains  an 
integer.  In  the  lower  left-hand  diagram,  the  parameter  k0  is  an  integer  while  the  cant  parameter  p0  is 
depicted  to  be  an  odd  multiple  of  1/2.  In  this  situation,  the  cant  can  never  be  fully  compensated  by  the 
integer  variable  v .  The  lower  right-hand  diagram  depicts  the  situation  where  both  k0  and  p0  are  odd 
multiples  of  1/2.  In  reality,  the  two  scale  offsets  can  fall  anywhere  between  ±  1/2.  In  Fig.  1  only  the 
extreme  offsets  are  depicted.  The  illustrations  demonstrate  the  functional  symmetry  derived  in  the  ear¬ 
lier  analysis. 

Hanning  Window 

To  study  the  effect  of  a  shading  window  function  on  the  CST,  the  popular  Hanning  window  is 
used.  In  this  case  the  window  function  wn  is  sin2  ( irn/N )  and  N/W  2.  Unfortunately,  a  closed- 
form  solution  is  not  available;  however,  by  using  Eq.  (10),  graphical  plots  of  the  characteristic  function 
were  computed  and  are  displayed  in  Fig.  2.  The  rationale  for  the  four  diagrams  is  identical  to  that 
described  for  the  case  of  the  rectangular  window  function. 

Cant  Selectivity 

A  study  of  the  characteristic  function  topology  for  the  two  window  functions  (Figs.  1  and  2) 
reveals  marked  differences  in  the  cant  selectivity.  In  the  case  of  the  rectangular  window  (commonly 
called  the  no  window  case  since  wn  —  1  over  the  integration  interval),  the  peak  of  the  contour  is  more 
sharply  defined.  The  variation  of  the  contour  over  the  x,y  plane  is,  however,  more  erratic.  At  the 
larger  values  of  y ,  the  spread  of  the  peak  along  the  x  axis  is  relatively  broad.  On  the  other  hand,  the 
topology  of  the  characteristic  function  for  the  Hanning  window  is  relatively  smooth,  and  it  displays  a 
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Fig.  1  —  Topology  of  the  CST  characteristic  function  Z{x,y)  for  a  rectangular  window.  The  upper  left-hand  diagram 
illustrates  the  case  when  k0  and  i/0  are  integers.  In  the  upper  right-hand  diagram,  i/0  is  an  integer  but  kQ  is  an  odd  multiple 
of  1/2;  that  is,  k0  occurs  at  a  bin  edge.  The  two  peaks  occur  for  y  =  v-  v0=  ±  1.  In  the  lower  left-hand  diagram,  k0  is  an 
integer  but  v0  is  an  odd  multiple  of  1/2;  that  is,  is  midway  between  two  cant  values  of  v.  The  two  central  peaks  are  equal 
and  occur  for  y  =  v-v0  =  ±1/2.  In  the  lower  right-hand  diagram,  both  k0  and  v0  are  odd  multiples  of  1/2. 
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Pig  2  —  Topology  of  the  CST  characteristic  function  Z(x^y)  for  a  Hanning  window.  The  rationale  and  peak  characteristics  of 
the  four  diagrams  are  the  same  as  described  in  Fig  1.  The  topology  for  the  shaded  window  is,  however,  more  smooth.  And, 
the  peak  of  the  topology  decays  less  rapidly  (along  the  line  of  minimum  phase  excursion  x=—y! 2)  as  the  absolute  value  of  y 
increases. 
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more  gradual  falloff  along  the  line  y  =—2x.  For  the  purpose  of  estimating  the  cant  of  a  linear 
trequency-slide  signal,  the  rectangular  window  would  prove  superior;  particularly,  in  a  noisy  environ¬ 
ment. 


To  obtain  a  quantitative  measure  of  the  cant  selectivity,  the  peak  value  of  the  characteristic  func¬ 
tion  over  x  was  computed  as  a  function  of  y  for  the  two  window  functions.  These  values  involve 
offsets  within  ±1/2  about  integer  values  of  x  and  y.  The  results  are  plotted  in  Fig.  3  as  a  function  of 
integer  values  of  the  cant  variable.  The  diagrams  display  the  expected  value  (small  circle)  and  the 
range  of  variation  of  Z(x,y )  due  to  the  random  offsets.  The  value  of  Z(x,y)  at  each  y  in  Fig.  3  are 
uniquely  dependent  on  the  specific  offsets.  Thus,  the  indicated  variations  along  the  abscissa  are  not 
independent  but  are  highly  correlated.  The  results  demonstrate  the  superiority  of  the  rectangular  win¬ 
dow  in  resolving  the  cant  of  a  narrowband  signal  in  a  noisy  background  environment. 
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Hg.  3  C  ant  selectivity  of  the  CST  characteristic  function  Z(x,y )  for  the  rectangular  and  Hanning  windows  The 
vertical  lines  depict  the  range  of  variation  in  CST  output  for  k0  and  uniformly  distributed  over  ±1/2  about  given 
integers.  The  specific  values  at  points  along  the  abscissa  are  not  independent  but  perfectly  correlated.  That  is,  given 
ihe  values  for  k0  and  v0 ,  the  CST  outputs  are  uniquely  determined  at  the  points  along  the  abscissa  and  will  fall 
within  the  range  indicated  in  the  figure. 


CST  STATISTICS  FOR  A  RANDOM  NOISE  SIGNAL 

Let  Fiik,v)  be  the  CST,  given  in  Eq.  (3),  of  a  sample  signal  randomly  selected  from  an  ensemble 
of  zero-mean  Gaussian  functions  with  standard  deviation  a.  Define  v  as  the  specific  value  of  „ 
whereby  F,(k,v)  <  Ft{k,i)  for  a  given  value  of  k.  Therefore,  F,(k,v)  and  Fj{k,i>)  are  samples  of  a 
random  function  whose  statistics  are  studied  in  the  following  paragraphs.  For  purposes  of  the  study, 
the  CST  size  is  assumed  to  be  very  large,  much  larger  than  the  range  over  which  the  cant  v  is  varied. 


Distribution  of  Ft  (A:,  v) 

For  a  Gaussian  signal,  it  is  known  that  Ft{k^o)  is  a  Rayleigh-distributed  random  variable  with 
mean  proportional  to  -TVila-  and  variance  equal  to  4/tt-I  times  the  square  of  the  mean.  Computer 
runs  over  a  large  set  of  Gaussian  inputs  have  verified  that  F,(k,v)  is  Rayleigh  distributed  and  identical 
for  each  k  and  v  (as  expected). 
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Distribution  of  Ft(k,v) 

To  study  the  distribution  of  the  peak  CST  for  a  Gaussian  noise  input  let 

Xj  =  Fj{k,vj),  (13a) 


where  vj  is  selected  from  the  set  {v}.  Consider  the  set  Xi>X2» . . .  ,x„  and  denote  \  =  Fj(k,i>)  as  a 
member  of  the  set;  such  that,  Xj  <  X  for  j  =  1,2 Now,  it  has  been  shown  that  when  the  distri¬ 
bution  functions  for  the  set  of  \j  are  identical  and  independent,  the  distribution  function  for  \  is  sim¬ 
ply  [9] 

F-  (z)  =  IFX  (z)  V  , 

where  F-(z)  is  the  probability  that  x  is  less  than  z  and  Fx(z)  is  the  probability  that  Xj  is  less  than  z  for 
a  given  vj.  Since  x  is  Rayleigh  distributed, 

Fk(z)  =  (1  -p)\ 

=  t(-l  )jnCjPj,  (14a) 

j- 0 

where  „C}  is  the  binomial  coefficient  and  p  is  the  complement  Rayleigh  distribution,  exp  (— z2/2cr2}, 
or  the  probability  that  z  <  x-  The  complement  distribution  function  for  y  as  a  function  of  p  is  there¬ 
fore 


.M 


=  np 


i  n  —  1  ,  (n  —  1)  (/i  —  2)  9 

1 - —  p  + - - - - pl  -  . . . 


2! 


3! 


(14b) 


Thus  when  p  is  less  than  0.2/(n-l),  the  complement  distribution  function  for  x  approximates  np. 

Curves  of  the  idealized  complement  distribution  function  are  shown  in  Fig.  4  for  several  values  of 
n .  Since  these  curves  are  based  on  independent  distribution  functions  for  the  set  Xj->  computer  runs 
were  made  to  determine  the  difference  between  the  actual  and  idealized  complement  distribution  func¬ 
tion  for  the  same  values  of  n .  This  proved  to  be  a  rather  formidable  task  due  to  the  large  sample  size 
required  to  achieve  accurate  measures  of  the  function.  However,  sufficient  data  were  achieved  to  show 
that  the  difference  is  quite  small.  In  every  case,  the  actual  probability  proved  to  be  slightly  less  than 
the  data  given  in  Fig.  4;  particularly  in  the  range  of  p  between  10-1  and  10~4.  The  results  for  n= 41 
(where  the  accuracy  is  higher  than  for  the  lower  values  of  n)  are  shown  in  Fig.  5.  The  results  for  lower 
values  of  n  are  comparable  to  those  shown  in  the  figure  for  values  of  p  less  than  about  10~2.  The 
difference  is  that  as  n  becomes  smaller,  the  curves  descend  to  about  the  same  minimum  values  for 
slightly  higher  values  of  p  before  merging  in  close  proximity  with  the  curves  shown.  As  a  conse¬ 
quence,  one  will  not  err  greatly  in  using  the  curves  shown  in  Fig.  4. 

Mean  and  Standard  Deviation  of  Fjik^v) 

The  probability  density  of  x  is  simply  the  derivative  of  the  distribution  function  (14a)  or 

fAz)  =  p(l-p)n-'  (15a) 

*  cr 

where  p  is  the  complement  Rayleigh  distribution  function.  The  various  moments  may  be  computed 
from  the  probability  density  function.  Carrying  out  the  indicated  operations  gives 

x/x  =  n  X,  n-iCj  ,  (15b) 
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PROBABILITY  p 


Fig.  4  -  Complement  distribution  function  of  the  CST  peak  output  for 
a  random  Gaussian  signal  as  a  function  of  the  complement  Rayleigh 
distribution  function,  p.  The  curves  assume  the  distribution  functions 
of  the  CST  at  each  u  are  independent.  At  small  values  oT  p  the  curves 
approximate  np ,  where  n  is  the  total  number  of  cant  values  in  the 
CST. 
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Fig.  5  —  Experimentally  determined  deviations  from  the 
complement  distribution  function  shown  in  Fig.  4  as  a  result  of 
dependence  between  the  cant  values  v\  n= 41.  The  deviation  for 
other  values  of 'n  closely  approximate  those  illustrated  for  p  less 
than  about  10  2.  As  n  gets  smaller,  however,  the  deviation 
tends  to  decay  and  reach  approximately  the  same  minimum  for 
slightly  larger  values  of  p .  It  may  be  concluded  that  one  will  not 
err  significantly  by  using  the  curves  shown  in  Fig.  4  for  the 
complement  distribution  function  of  the  CST  peak  output. 
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and 


°S/<rx  = 


4  n 


I 


J- o 


(-1V 

o+o2 


n-lCj  —n-(x/x;)2 


1/2 


(15c) 


where  \  and  are  the  mean  and  standard  deviation  of  the  Rayleigh  distributed  CST  statistic 
Xj  =  Fj  (k,vj). 

Figure  6  is  a  plot  of  the  above  functions  along  with  experimental  measures  obtained  from  the 
CST  by  using  a  Gaussian  noise  signal.  The  experimental  data  reflect  the  dependence  or  covariance 
between  the  Xj  sample  increments. 


NUMBER  OF  CONTIGUOUS  CANT  VALUES  IN  CST,  n 

Fig.  6  —  Mean  and  standard  deviation  of  the  peak  CST 
output  for  a  random  Gaussian  signal  as  a  function  of  n , 
relative  to  the  values  for  n—1.  The  theoretical  curves 
assume  independence  between  the  cant  values,  while  the 
dashed  curves  were  obtained  experimentally  by  averaging 
about  250,000  samples  per  point.  Deviation  from  the 
theoretical  curves  is  attributable  to  the  dependence 
between  cant  values. 


Characteristics  of  the  Variable  v 

Experimental  measures  were  made  to  determine  the  distribution  of  the  statistic  v  for  a  random 
Gaussian  signal.  This  variable  is  the  member  of  the  set  M  that  maximizes  the  CST  for  a  given  input 
signal  sample.  The  observed  distribution  of  v  is  uniform  over  the  range  of  cants  except  at  the  two  end 
points  where  the  probability  is  higher  than  the  mean.  This  implies  that  the  covariance  between  the  cant 
parameters  are  nonzero.  In  this  event,  maximum  CST  values,  which  would  occur  for  v  values  outside 
of  the  range  of  cants  employed  in  the  CST,  will  increase  the  probability  that  the  maximum  will  occur  at 
the  nearest  end  point. 
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To  verify  the  above  hypothesis,  the  normalized  covariance  between  v.  and  y, ,  from  Ea  (13a) 
was  computed  experimentally.  Figure  7  shows  the  results  for  both  the  rectangular"  Hanning 

dow  °and  fn  ear  y’  3  Slgmflcant  correlat*on  exists  between  adjacent  values  of  v  for  the  rectangular  win- 
, fgreater  separation  of  v  values  in  the  case  of  the  Hanning  window.  This  further  veri- 
use  o  o  her  than  the  rectangular  window  is  disadvantageous  in  the  case  of  the  CST  Com 
Putmg  ^  value  and  standard  deviation  of  the  CST,  using  a  sequence  of  cant  values  separated  4 
or  greater  apart,  resulted  in  values  that  approached  the  theoretical  curve  shown  in  Fig.  6. 


Fig-  7  Normalized  covariance  between  cant  values  as  a 
function  of  the  cant  separation  A*/  for  rectangular  and 
Hanning  windows. 
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Experimental  Parameters 

The  desired  distributions  were  measured  experimentally  by  using  Eq.  (16)  in  Eq,  (3),  for  both  a 
rectangular  and  a  Hanning  window.  The  CST  parameters  were;  N  =  1024,  k  =  200,  and  v  —  —  20, 

•  •  • ,  0,  •  •  • ,  20.  The  signal  parameters  included:  k0  =  200  and  200.5,  v0  =  —0.5,  0,  9.5,  and  10;  and 
inband  signal-to-noise  ratio  r  =  20,  10,  6,  3,  0,  -3,  and  -6  dB.  Using  Gaussian  noise  statistics,  more 
than  20,000  samples  of  the  CST  statistic  were  obtained  for  each  set  of  parameters  in  measuring  the  v 
distributions.  The  v  distributions  for  Gaussian  noise  only  were  determined  earlier. 

Distributions  of  the  Peak  Cant  v 

As  expected,  the  distributions  of  v  are  independent  of  the  parameter  v0  as  long  as  v0  is  an  interior 
cant  within  the  range  of  the  cants  vJt  and  is  not  close  to  an  end  value.  The  distributions  are  dependent 
on  the  offsets  of  both  k0  and  v0  from  their  nearest  integer  values.  Figure  8  shows  the  distributions  of  v 
as  a  function  of  the  inband  signal-to-noise  ratio  r,  for  a  rectangular  window  and  for  fractional  offsets  of 
the  parameters  k0  and  v0.  The  upper  left-hand  diagram  shows  the  distributions  when  both  k0  and  v0 
are  integers.  The  upper  right-hand  diagram  illustrates  the  effect  of  Atq  being  offset  by  1/2;  that  is  to  the 
edge  of  the  k  bin.  The  lower  left-hand  diagram  illustrates  the  case  where  v0  is  midway  between  two 
integer  cants.  And  the  lower  right-hand  diagram  illustrates  the  case  where  both  ko  and  are  offset  by 
1/2.  From  the  figure,  one  can  readily  see  that  a  cant  shift  of  one  nearly  compensates  for  the  degrada¬ 
tion  induced  by  the  signal  frequency  occurring  at  the  edge  of  the  spectral  bin.  Further,  the  degradation 
resulting  from  the  signal  cant  being  midway  between  two  integer  cant  values  is  partially  compensated 
when  the  signal  frequency  is  offset  from  the  bin  center.  In  practice,  the  values  of  /cq  and  v$  can  be 
expected  to  be  uniformly  distributed  over  ±1/2  from  an  integer  value. 

Figure  9  displays  the  distributions  of  v  as  a  function  of  the  inband  signal-to-noise  ratio  r,  for  a 
Hanning  window  and  for  fractional  offsets  of  the  two  signal  parameters.  The  rationale  of  the  four 
diagrams  is  identical  to  that  given  for  the  case  of  the  rectangular  window.  The  distributions  in  the  k 
bin  do  not  differ  greatly  from  those  realized  with  the  rectangular  window.  The  peak  values  are  approxi¬ 
mately  the  same,  although  the  distributions  about  the  peak  values  are  somewhat  broader.  This 
broadening  effect  about  the  peaks  can  be  attributed  to  the  broader  bandwidth  of  the  CST  for  a  Hanning 
window.  As  evidenced  from  Figs.  1  and  2,  the  v  distributions  for  values  of  k  remote  from  A:0  can  be 
significantly  different  for  the  two  window  functions. 

Cumulative  Distributions  of  the  Peak  Cant  v 

The  cumulative  distribution  of  the  peak  cant  is  defined  as  the  sum  of  the  probabilities  at  0  and 
about  either  side  of  zero  on  the  abscissa  scale  in  the  diagrams  in  Figs.  8  and  9.  That  is,  the  cumulative 
distribution  at  a  scale  value  of  m  is  simply  the  probability  at  the  scale  value  0  plus  the  sum  of  the 
probabilities  at  both  positive  and  negative  scale  values  up  to  and  including  m.  Using  this  definition,  the 
cumulative  distributions  of  the  peak  cant  for  the  rectangular  window  are  shown  in  Fig.  10.  The 
rationale  of  the  four  diagrams  is  the  same  as  that  described  for  the  earlier  distributions.  In  the  upper 
left-hand  diagram,  the  cumulative  distribution  for  r  equal  to  20  dB  coincides  with  the  vertical  ordinate 
axis,  and  consequently  is  not  perceived.  It  is,  however,  visible  on  the  remaining  three  diagrams  as 
illustrated  in  the  upper  right-hand  diagram.  As  may  be  seen,  the  probability  rises  rapidly  at  a  scale  of 
one,  after  which  it  rises-  less  rapidly  with  increased  values  of  the  abscissa  scale.  The  cumulative  distri¬ 
bution  for  a  noise-only  signal  ( r=-°°dB )  is  simply  a  straight  line  as  expected.  These  curves,  as  well  as 
the  earlier  distributions  represent  the  results  for  41  contiguous  values  of  the  cant  variable. 

The  cumulative  distributions  of  the  peak  cant  for  the  Hanning  window  are  displayed  in  Fig.  11, 
with  the  same  format  as  employed  in  Fig.  10.  The  curves  for  the  Hanning  window  closely  approximate 
those  for  the  rectangular  window.  They  generally  rise  more  steeply  at  an  abscissa  value  of  one;  after 
which  they  ascend  more  gradually.  This  behavior  is  attributed  to  the  broader  peaks  of  their  distribu¬ 
tions. 
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V  -  VO  -  1/2  v  -  „o  -  1/2 

Fig.  8  —  Distributions  of  the  peak  cant  v  for  a  ramp-frequency  signal  in  random  noise,  a  rectangular  window,  and 
selected  values  of  inband  signal-to-noise  ratio  r.  The  upper  left-hand  diagram  illustrates  the  case  when  the  signal 
parameters  k0  and  p0  are  integers.  In  the  upper  right-hand  diagram,  v0  is  an  integer  but  kQ  is  an  odd  multiple  of  1/2; 
that  is,  k0  is  at  the  bin  edge.  In  the  lower  left-hand  diagram,  k0  is  in  the  center  of  the  k  bin  but  is  an  odd  multiple 
of  1/2;  that  is,  u0  is  midway  between  two  cant  values.  In  the  lower  right-hand  diagram,  both  k0  and  u0  are  odd 
multiples  of  1/2.  The  number  of  contiguous  cant  values,  in  all  cases,  is  41. 
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Fig.  9  Distributions  of  the  peak  cant  i  for  a  ramp-frequency  signal  in  random  noise,  a  Hanning  window,  and 
selected  values  of  inband  signal-to-noise  ratio  r.  The  rationale  for  the  four  diagrams  is  the  same  as  described  in  Fig. 
8.  The  number  of  contiguous  cant  values,  in  all  cases,  is  41. 
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|  v  —  vq  —  1/2 1  \v  —  vo  —  1/2 1 


Fig.  10  —  Cumulative  distributions  of  the  peak  cant  v  for  a  ramp-frequency  signal  in  random  noise,  a  rectangular 
window,  and  selected  values  of  inband  signal-to-noise  ratio  r.  The  rationale  for  the  four  diagrams  is  the  same  as 
described  in  Fig.  8.  The  number  of  contiguous  cant  values,  in  all  cases,  is  41  The  cumulative  distributions  are 
computed  from  the  data  in  Fig.  8  by  summing  the  probabilities  for  the  abscissa  scale  value  of  0  and  the  absolute  scale 
values  about  0  up  to  the  indicated  scale  value. 
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I  v  ~  V ol  I*  ~  vq\ 


v  —  vq  —  '\/2\  |  ^  —  i^o  —  1/2 1 


Fig-  11  —  Cumulative  distributions  of  the  peak  cant  v  for  a  ramp-frequency  signal  in  random  noise,  a  Hanning 
window,  and  selected  values  of  inband  signal-to-noise  ratio  r.  The  rationale  for  the  four  diagrams  is  the  same  as 
described  in  Fig.  8.  The  number  of  contiguous  cant  values,  in  all  cases,  is  41,  The  cumulative  distributions  are 
computed  from  the  data  in  Fig.  9  by  summing  the  probabilities  for  the  abscissa  scale  value  of  0  and  the  absolute  scale 
values  about  0  up  to  the  indicated  scale  value. 
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It  is  well  to  remember  that  all  of  the  distributions  shown  in  Figs.  8  through  1 1  are  truly  discrete 
values  and  not  continuous.  Straight  lines  between  the  discrete  values  of  probability  have  been  drawn 
for  convenience  in  displaying  the  results. 

SUMMARY  AND  CONCLUSIONS 

The  conventional  sectionalized  Fourier  transform  (SFT)  accumulates  signal  energy  in  narrow 
spectral  channels  or  frequency  bins  that  are  constant  over  the  integration  time  of  the  transform  It  has 
been  shown  that  by  extending  the  phase  kernel  of  the  SFT,  signal  energy  over  narrow  spectral  channels 
that  vary  dynamically  over  the  integration  period  can  be  accumulated  to  achieve  optimum  detection 
sensitivity  for  a  narrowband  signal  whose  instantaneous-frequency  dynamics  are  known.  When  the 
spectral  dynamics  of  the  signal  are  unknown,  the  parameters  of  the  generalized  kernel  may  be  varied  to 
obtain  estimates  of  the  input  signal-frequency  dynamics.  The  canted  spectral  transform  (CST)  which 
employs  a  quadratic  phase  kernel  (or  Fourier-Fresnel  kernel),  is  an  initial  step  toward  achieving  a  gen¬ 
eralized  spectral  transform. 

The  statistical  properties  of  the  CST  reveal  that  shaded  window  functions  are  generally  inferior  to 
the  rectangular  window  (commonly  referred  to  as  no  window)  in  discriminating  the  cant  of  signals 
whose  instantaneous  frequency  varies  linearly  with  time.  Although  the  primary  utility  of  the  CST  is  in 
detecting  signals  with  linearly  varying  frequency  characteristics,  it  also  has  merit  as  a  conventional  spec¬ 
trum  analyzer  by  appropriately  exploiting  the  distribution  characteristics  of  the  peak  cant  variable.  The 
statistical  properties  of  the  CST  for  both  signal  and  noise  provide  a  basis  for  determining  the  perfor¬ 
mance  of  the  transform  in  practical  applications. 
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